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ABSTRACT 

We measure the cross-power spectrum between galaxy density from Canada-France- 
Hawaii-Telescope Lensing Survey (CFHTLenS) catalogues and gravitational lensing 
convergence from Planck data release 1 (2013) and 2 (2015). We investigate three 
main galaxy samples: 18.0 < tAB < 22.0, 18.0 < tAB < 23.0, 18.0 < tAB < 24.0 
in the redshift range 0.2 < z < 1.3 in each of the four CFHTLenS wide fields. 
By comparing the measured cross-spectrum with model predictions, linear galaxy- 
dark matter biases of 6 = 0.82to;23J 0-83lo;i[g, 0.82^Q;i[4 are inferred at significances of 
3.5,4.5, 5.6tT using the Planck 2015 release. These measurements are marginally consis¬ 
tent with biases derived from galaxy-galaxy auto-correlations: h = 1.15to oi, l.OSlgY 
and 0.96;)!Y respectively. Using the 2013 Planck release, we obtain biases of 6 = 
1.33lo;2g, 1.19 ^q; 23> showing significant differences between the releases. 
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1 INTRODUCTION 

An approach to measuring the distribution of dark mat¬ 
ter in our Universe is by using weak gravitational lensing, 
where light is deflected by mass between the source and 
the observer. In the case of cosmic microwave background 
(CMB) weak lensing, the paths of photons emitted from 
the last scattering surface become perturbed by the large- 
scale structure in the Universe amidst their journey to the 
observer. This changes the statistics of hot/cold spots size 
distribution, leading to a change in the local angular power 
spectrum on the scale of the primary acoustic peaks. ( |Lewis| 
I&: Challinor| 20 d 6 1 . Since the degree of deflection is propor¬ 
tional to the integrated potential that the photon traverses, 
weak lensing is a sensitive probe for studying the large-scale 
structure between the CMB and us. 

The clustering properties and statistics of dark mat¬ 
ter halos have been extensively studied through simula- 


20101. The ratio between the clustering amplitude of dark 


tions I Boylan-Kolchin et al.|2d09 Klypin et al.|20lT |. Since 
these dark matter halos form in peaks in the peaks of the 
background density distribution, they are clustered more 


strongly than the background matter distribution (Kaiser 


1984 Bardeen et al. 19861. Although galaxies live inside 


dark matter halos, the formation processes and their dynam¬ 
ics make the distribution of galaxies non-local and possibly 
stochastic relative to the halo center, making the cluster¬ 


ing characteristics different (Mo & van den Bosch & White 


matter and galaxies is called linear galaxy bias. 

Here we cross-correlate the galaxy number density 
from the Canada-France-Hawaii Telescope Lensing Survey 
(CFHTLenS) ||Heymans et al.|[M^ |Erben et al.][2^ ^ 


[lander et al.||2014| ) and the all-sky convergence map from 
Planck (Planck XV 2015, hereafter PXV2015) to measure 
galaxy biases. Similar studies include Planck lensing cross- 
correlated with NRAO VLA Sky Survey (NVSS) galax¬ 
ies, MaxBCG clusters, SDSS LRGs and WISE (Planck 
XVH 2013, hereafter PXVH2013), Wilkinson Microwave 
Anisotropy Probe (WMAP) lensing cross-correlated with 
NVSS galaxies ( [Smith, Zahn &: Dore|2007| |, WMAP lensing 
cross-correlated with LRGs and quasar from SDSS (Hirata 
et al.|2008 1 Blanco Cosmology Survey (BCS) galaxies cross- 


correlated with South Pole Telescope (SPT) lensing (|Bleem[ 


et al. 20121, Sloan Digital Sky Survey (SDSS) quasar 


maps cross-correlated with Atacama Cosmology Telescope 
(ACT) lensing ( [Sherwin et al.[[2012[ ). Wide-field Infrared 
Survey Explorer (WISE) selected quasars cross-correlated 
with SPT lensing ( Geach et al.|201^ , Herschel/SPIRE cos¬ 
mic infrared background (CIB) cross-correlated with SPT 
lensing (Holder et al. |2013| , SPTpol E-mode polariza¬ 
tion cross-correlated with estimates of the lensing poten¬ 
tial from Herschel/SPIRE CIB ( [Hanson et al.[[2013[ ), Po- 
larbear experiment CMB polarization cross-correlated with 
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Herschel/SPIRE CIB (Polarbear collaboration 20141 and re¬ 
constructed convergence from ACTpol cross-correlated with 
CIB measurements from Planck (van Engelen et al.|2014b|. 
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2 Y. Omori and G. Holder 


In our study, we present both the cross-correlations be¬ 
tween galaxy overdensity and lensing, and the galaxy auto¬ 
correlations using the same galaxy data, and compare the 
linear galaxy biases that we obtain. 

In this paper we assume flat ACDM cosmology {yim = 
0.28, Ha = 0.72, Ho = lOOh km s“^Mpc“^ and erg = 
0.82, Us = 0.96) with h = 0.70. All magnitudes are given 
in the AB system. 


2 THEORY 


CMB lensing convergence is calculated by integrating the 
matter fluctuation in the line-of-sight direction: 

K{n) = J dxW'"{x)S[xn,z{x)], ( 1 ) 

where n is the line-of-sight direction, x is the line-of-sight 
comoving distance, <5 is the fractional dark matter density 
fluctuation and is the lensing kernel: 


2 V c / a{x) XCMB 


( 2 ) 


which describes the efficiency of lensing at a given redshift 
when multiplied by dx/dz. 

Under the assumption of the galaxy distribution trac¬ 
ing out the peaks of dark matter density fluctuations, the 
overdensity of galaxies can be described as: 


fl(n) = / dxW^{xY[xn,z(x)]- 


(3) 


The distance kernel iy^( |Sherwin et al.|2012 l: 

1 dz dN{z 


W^ix) = 


J dz 


, dN{z') 
dz' 


dx dz 


b{x) + MB{x) (4) 


describes the distribution of galaxies as a function of redshift 
in terms of co-moving distances. The second term in equa- 
tion[^is the magnification bias, where the observed number 
density of galaxies are altered due to lensing caused by mass 


between the foreground galaxy and the observer (Turner & 
Qstriker &: Gott|[l984 \ is given by: 

2 




Ho 

c 


(l + 2(x))5(x) (5s-2), 

(5) 


with 


and 


/■XCMB -y' _ 

gix) = X dx'^ 

Jv X 


x' -X dN/dz dz 
'jdz'dN/dz'dx'’ 


dlogioiV(< m) 
dm 


( 6 ) 


(7) 


for a galaxy sample with limiting magnitude mum (Villum- 
sen||1995 1. As shown in figure]^ the amplitudes of the magni¬ 
fication biases for each sample are significantly smaller than 
the dN/dz term. 

The cross-spectrum between convergence and galaxy 
overdensity under the Limber approximation is the power 


spectrum weighted by the lensing and galaxy kernels I Lim- 
ber||1953 Kaiser|[l992 1: 


= / dz^^WYx}W^ix)p(k=^,zy 


( 8 ) 



z 


Figure 1. Light green, Turquoise, Dark green: Redshift 

distribution dn/dz for the 18.0 < 2 ab < 22.0, 18.0 < ^ab < 23.0, 
18.0 < 2AB < 24.0 samples in the range 0.2 < 2 < 1.3. The 
distribution is normalized such that the integral f dz(dNfdz) = 1. 
Dashed line: Convergence kernel multiplied with (dxjdz) 
and arbitrary normalization applied for visualization purposes to 
demonstrate the effectiveness of the lenses at various redshifts. 
Lower panel: Amplitude of magnification bias for each of the 
galaxy samples. 


where P{k,z) is the non-linear matter power spectrum at 
redshift 2 :. As a comparison, we also compute the galaxy 
overdensity auto-spectrum given by: 


/' 


( 7 ^= / dz'^^W^{xfp{k=^,z 


(9) 


The model spectra in equations are computed using 
the non-linear power spectrum from CAMB ([Lewis fc Bridle 
2002 1 and revised Halofit (Smith et al. 2003[ Takahashi 
et al.|20T2 1. The redshift distributions dN/dz for each sub¬ 
sample are produced by averaging the individual photomet¬ 
ric redshift probability distributions for each galaxy rather 
than simply using the best-fit estimates for each galaxy. 


3 GALAXY MAPS 

In this study, we use the galaxies from the CFHTLen^ 
galaxy survey. CFHTLenS is part of the wide component 
of the Canada-France-Hawaii Telescope Legacy Survey 
(CFHTLS), which consists of 4 fields centered at 2:18:00 
-07:00:00, 08:57:49 -03:19:00, 14:17:54 -b54:30:31, 22:13:18 
-1-01:19:00, each ranging from 23 — 64 deg^ with a total of 

^ http://cfhtlens.org/ 

^ http://www.cfht.hawaii.edu/Science/CFHTLS/ 
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154 deg^ and full multi-colour depth of zab = 24.7 ( Heyman^ 


et al.|20l2 Erben et al.|2013 Velander et al.|2014 1. We limit 

our analysis to galaxies in the redshift range 0.2 < 2 < 1.3, 
which is confirmed by Heymans et al. (20121 to have a 


photometric redshift distribution that resembles the mea¬ 
sured spectroscopic redshift distribution. Galaxies selected 
with lAB < 24.5 in this redshift slice have a scatter of 
0.03 < (TAz < 0.06 (where ctaz is the variance in the dif¬ 
ference between the photometric and spectroscopic redshifts 
{zp — 2 s )/(1 + 2 s )) with 10% of the galaxies classified as out¬ 
liers ( Hildebrandt et al.|2007 Benjamin et al.|2013 1. Within 
this redshift slice, three samples with 18.0 < iab < 22.0, 
18.0 < iAB < 23.0, 18.0 < lAB < 24.0 are extracted, re¬ 
sulting in galaxy catalogues with mean number densities 
n = 3.4, 7.7,15.4 per square arcminute respectively. The 
lower magnitude cut is applied to reduce dispersion and 
bias due to photometric redshift estimates in addition to 
removal of possible stars that passed through initial reduc¬ 
tion pipeline classifications (star_flag = 0 and mask = 0, 


descriptions of each flag can be found in Erben et al. (20131). 


Using these catalogues, the fractional galaxy overden¬ 
sity maps at arcminute resolution are produced by: 


5ij — 


Nij — {N)wi^ 
{N)wij 


( 10 ) 


for each i-th and j'-th pixel, where Wij is the pixel-by-pixel 
weight map produced by degrading the arcsecond resolution 
masks to match the same resolution and 




( 11 ) 


4 CONVERGENCE MAPS 

We use the observed and simulated convergence maps from 
the Planck 2015 releas^ (PXV2015) and compare the re¬ 
sults with the 2013 releas^ (PXVII2013). The 2015 con¬ 
vergence map is produced by applying a quadratic estima¬ 
tor ( [Qkamoto Hu]|2003[ ) on a SMICA component sepa¬ 
rated map, which is synthesized by cleaning the foregrounds 
and combining all nine frequency bands. The Galaxy and 
point sources are masked in the process and the bands are 
restricted in the range of 100 < £ < 2048 (PXV2015). The 
2013 map is released as a. f map is generated by combin¬ 
ing the 143 and 217 GHz channels with dust distribution 
subtracted using the 857 GHz channel. The map also masks 
point sources included in the ERCSC, SZ and PCGS catalogs 
(PXVH2013). The map is transformed into a convergence k 
map by taking the transform: 

(12) 

where TZf is a normalization factor. We combine the Planck 
2013 and 2015 masks and apply this to both convergence 
maps for consistency. Additional masking is not applied be¬ 
yond the masks released by the Planck collaboration. Field 

^ http://irsa.ipac.caltech.edu/data/Planck/release_2/all-sky- 

maps/maps/component-maps/lensing/COM_CompMap_Lensing_ 

2048_R2.00.tar 

http://irsa.ipac.caltech.edu/data/Planck/release_l/all-sky- 
maps/previews/COM_CompMap _Lensing_2048-Rl. 10/index.html 
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Figure 2. Cross-correlation of 18.0 < jab < 24.0 sample for 
W1,W2,W3,W4, and the average of the four fields with the best 
fit theory model. 


regions were extracted by multiplying this all-sky mask with 
galaxy catalogue masks. 

Biases in the lensing reconstruction due to the ther¬ 
mal Sunyaev-Zeldovich (tSZ) effect or cosmic infra-red back¬ 
ground (GIB) fluctuations are expected to be small with 
Planck’s angular resolution (van Engelen et al. 2014a Os- 


borne, Hanson fc Dore|[Ml4K 


5 CROSS AND AUTO-SPECTRUM 


The cross-power spectrum is calculated using the flat-sky 
approximation and multiplying the maps in Fourier space: 

= fiReiim - ( 13 ) 


where / is the normalization factor N 1 N 2 / and 1C, 

Q are the Fourier transforms of the convergence and the 
galaxy overdensity maps multiplied with the weight map. 
The auto-spectrum is calculated in a similar way but re¬ 
placing Q with 1C and subtracting the shot noise contribu¬ 
tion = l/(n), where (n) is the average number density 

of galaxies per steradian. 

Uncertainties are calculated by cross-correlating the 
galaxy maps with 100 simulated Planck convergence maps 
(released by the Planck collaboration) reduced through the 
same pipeline, and calculating the variance for each bin. 
Using these uncertainties, the weighted averages for each 
individual bin are calculated. 

Since the uncertainties in the galaxy auto-spectrum 
cannot be derived in the same way, we resort to an inter¬ 
nal block-jackknife method often used in 2-point correlation 
function analysis, where we divide each of the fields into 
64 sub-regions and the jackknife samples are then created 
by omitting each of these sub-samples ( [Zehavi et al.|[^05| 
Norberg et al.|2009 1. The covariance matrix is estimated by: 
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^ ■ {C,)UC7 - {C,))„ (14) 

n=l 

while the covariance matrix for the cross-correlation is cal¬ 
culated by: 

■^maps 

s”" = - E - (C^MCr - {Ci))j. (15) 

4’maps 

^ m = l 

Using these covariance matricies, the best ht biases are de¬ 
rived by minimizing the x^'- 

- Cf)j, (16) 

ij 

where X = ng or gg and i,j are the bin numbers. The 
signihcance of detection is defined by Ax^ = X^i^ ~ 0) ~ 
xLinib)- 


6 RESULTS 

The cross-spectra for each of the helds using the 18.0 < 
lAB < 24.0 sample are shown in hgure and the different 
magnitude samples are shown in figure I^The best ht values 
b and their associated x^ values are listed in table where 
the x^ values were calculated using 20 bins in the range of 
^ = 50 - 1900. 

The best ht biases for the Planck 2015 release are found 
to be 6 = 0.82l);;23,0.834o;l8,0.82lg;l^ at signihcances of 
3.5, 4.5, 5.6 ct respectively. The minimum x^ for the hts are 
25.4, 18.6, 15.0 with a probability-to-exceed (PTE) of 0.15, 
0.48, and 0.72. These best ht biases are signihcantly lower 
than the biases found from the Planck 2013 release: h = 
1.33to'28> 1'194 o'23i l-16to'ig. Similarly, the best ht biases 
from C®® are found to be l.lhto'oi) 1.084o;oii 0'964o'.ol, with 
X^ of 14.1, 22.1, 38.0. Sample correlation matricies are shown 
in hgure 

These results suggest that unlike quasars, which are 
strongly biased tracers of dark matter (see Geach et al. 2012, 
Sherwin et al. 2012), optically-selected galaxies in these 
magnitude ranges on average are unbiased tracers of dark 
matter huctuations. 


stronger constraints on the value of the bias, these cross¬ 
correlations provide a more direct measurement of the re¬ 
lationship between galaxies and two-dimensional projected 
mass, because the degree of photon path dehection is di¬ 
rectly correlated with the depths of foreground gravita¬ 
tional potentials that it passes through. Furthermore, cross¬ 
correlations are less sensitive to the complexity of sur¬ 
vey masks and any hidden systematic errors that could be 
present in any single survey. 

The power of cross-correlations is yet to be used at full 
potential. With the advent of upcoming wide galaxy sur¬ 
veys probing to fainter magnitudes, the sky coverage will 
improve the statistics while the depth will provide sufficient 
galaxy number density to maintain signal above shot-noise. 
For example, similar cross-correlation analyses will be per¬ 
formed between the Dark Energy Survey (DES) galaxies - 
SPT/SPTpol lensing and Hyper Supreme Camera (HSC) 
survey galaxies - ACT/ACTpol lensing, which will have sev¬ 
eral thousand square degrees of overlap. With such large 
areas and fainter magnitude limits, the signal to noise is 
expected to be an order of magnitude larger. 
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7 CONCLUSION 

We have presented the results of cross-correlations between 
galaxy number density from CFHTLenS and CMB lensing 
convergence from Planck, along with galaxy-galaxy auto¬ 
correlations. We obtain the dark matter-galaxy bias using 
these two independent methods. For the cross-correlations 
of the 3 magnitude samples, we obtain best fit biases of 6 = 
0.82+g;2t0.834o;ls,0.82t^;;;| in the range of 50 < ^ < 1900 
using the 2015 release and b = l.OOj’^Q'Jg, 1-191 q; 23> 
using the 2013 release. From the galaxy auto-correlations, 
we obtain biases of 6 = 1.154o'oii l-OSto'oii 0-964g;oii which 
settle between the values obtained from the two data releases 
and agrees well with previously established results that indi¬ 
cate linear galaxy bias of ~ 1 at large scales (see for example 
Gaztanaga & Frieman 1994). 

Although the galaxy-galaxy auto-correlations place 
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18.0 < ^ab < 22.0 18.0 < ^ab < 23.0 18.0 < ^ab < 24.0 
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Figure 3. Cross-spectrum for the 18.0 < 2 ab < 22.0, 18.0 < 2ab < 23.0, 18.0 < ^ab < 24.0 galaxy samples cross-correlated with 
Planck convergence (filled bars: 2015, empty bars: 2013) (upper row) and the auto-spectrum using the same galaxy samples (lower row). 
The points with errorbars shown are the coarsely binned calculations for displaying purposes, and the the fits to the 20 bins shown 

in the background. 
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Figure 4. Correlation matrices between bins for cross-correlation and auto-correlation. 
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Sample 

N 

(n) [arcmin 

fe (Cf ) 


^2015 (C”®) 

X^ 

b20i3 (C;«) 

x^ 

18 < zab < 22 

1444906 

3.4 

1 1 c:+0.02 

-*-•-*^^-0.01 

14.1 

0-821^23 

25.4 

1 qq+0.29 

15.4 

18 < 2ab < 23 

3268383 

7.7 

1 r>o + 0.01 

22.1 

0.83«;« 

18.6 

1 IQ+0-23 

15.9 

18 < zab < 24 

6529185 

15.4 

n Qfi+0-01 

u-yo_o 01 

38.0 

0.82l°;i« 

15.0 

^•-*-'^-0.18 

15.4 


Table 1. Table of best fit linear bias and values calculated using 20 bins between 50 < € < 1900. 
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